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Abstract. Landen transformations are maps on the coefficients of an integral 
that preserve its value. We present a brief survey of their appearance in the 
literature. 



To Henry, who provides inspiration, taste and friendship 
1. IN THE BEGINNING THERE WAS GAUSS 

In the year 1985, one of us had the luxury of attending a graduate course on 
Elliptic Functions given by Henry McKean at the Courant Institute. Among the 
many beautiful results he described in his unique style, there was a calculation of 
Gauss: take two positive real numbers a and b, with a > b, and form a new pair 
by replacing a with the arithmetic mean (a + 6)/2 and b with the geometric mean 
\fab. Then iterate: 

a n + b n 



(1-1) a n+ i = n - , b n+1 = \J a n b n 

starting with clq = a and bg = b. Gauss [29] was interested in the initial conditions 
a = 1 and b = The iteration generates a sequence of algebraic numbers which 
rapidly become impossible to describe explicitly, for instance, 



(1.2) a 3 = — ^(1 + ^2) 2 + 2V2^2 ^1 + ^2 

is a root of the polynomial 

G(a) = 16777216a 8 - 16777216a 7 + 5242880a 6 - 10747904a 5 
+942080a 4 - 1896448a 3 + 4436a 2 - 59840a + 1. 

The numerical behavior is surprising; a^ and be agree to 87 digits. It is simple to 
check that 

(1.3) lim a„ = lim b n . 

n — >oc n — >oo 

See (|6.1I) for details. This common limit is called the arithmetic-geometric mean 
and is denoted by AGM(a, b). It is the explicit dependence on the initial condition 
that is hard to discover. 

Gauss computed some numerical values and observed that 

(1.4) an - b n - 1.198140235 
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and then he recognized the reciprocal of this number as a numerical approximation 
to the elliptic integral 

(1.5) / - " 



It is unclear to the authors how Gauss recognized this number - he simply knew 
it. (Stirling's tables may have been a help; |12 contains a reproduction of the 
original notes and comments.) He was particularly interested in the evaluation of 
this definite integral as it provides the length of a lemniscate. In his diary Gauss 
remarked, 'This will surely open up a whole new field of analysis' [271 115j . 

Gauss' procedure to find an analytic expression for AGM(a, b) began with the 
elementary observation 

(1.6) AGM(a, b) = AGM ' " 



and the homogeneity condition 

(1.7) AGM(Aa, Xb) = A AGM (a, b) . 

He used ([OIJ) with a = (1 + \fk) 2 and b = (1 — \fk) 2 , with < k < 1, to produce 

AGM(1 + k + 2Vk, 1 + k- 2Vk) = AGM(1 + fc, 1 - k). 

He then used the homogeneity of AGM to write 

AGM(1 + k + 2\/fc, 1 + k - 2\fk) = AGM((1 + ft)(l + fc*), (1 + fc)(l — k*)) 

= (1 + fc)AGM(l + k*,l- k*), 

with 

This resulted in the functional equation 

(1.9) AGM(1 + jfe, 1 - fc) = (1 + fc) AGM(1 + fc*, 1 - fc*). 

In his analysis of (|1.9p . Gauss substituted the power series 



i 00 

into (|1.9p and, solved an infinite system of nonlinear equations, to produce 



(1.11) a n = 2 



-2n 



2n\ 2 



Then he recognized the series as that of an elliptic integral to obtain 

1 2 r /2 dx 

(1.12) 



AGM(1 + fc, 1 - fc) 7T J y/i-k? sm 2 x' 
This is a remarkable tour de force. 

The function 

dx 



(1.13) 



K(k) = f 
Jo 



yl — fc 2 sin 2 x 
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is the elliptic integral of the first kind. It can also be written in the algebraic form 

»i i, 

(1.14) K(k) 



io V(l-^)(l-fc 2 ^ 2 ) 
In this notation, (|1.9|) becomes 
(1.15) K(k*) = (l + k)K{k). 

This is the Landen transformation for the complete elliptic integral. John Landen 
[37j . the namesake of the transformation, studied related integrals: for example, 

i dx 

(1.16) 

He derived idcntitcs such as 



o v/:r 2 (l-z 2 ) 



(1.17) k = e\Je 2 - 7T , where e := / y/l - sin 2 9 d6 



r/2 



proven mainly by suitable changes of varibles in the integral for e. In 52J the reader 
will find a historical account of Landen's work, including the above identities. 

The reader will find in [T5] and [12] proofs in a variety of styles. In trigonometric 
form, the Landen transformation states that 

^ t2 dO 

(1.18) G(a,b)= 1 



\J o? cos 2 9 + b 2 sin 2 9 

is invariant under the change of parameters (a, b) i— > M^^jVa&J. D. J. Newman 
[46] presents a very clever proof: the change of variables x = b tan 6* yields 



(1.19) G{a,b)= l -J" 



v/(a 2 +x 2 )(6 2 +x 2 ) 



Now let x i— > a; + Vx 2 + ab to complete the proof. Many of the above identities can 
now be searched for and proven on a computer |12j . 

2. An interlude: the quartic integral 

The evaluation of definite integrals of rational functions is one of the standard 
topics in Integral Calculus. Motivated by the lack of success of symbolic languages, 
we began a systematic study of these integrals. A posteriori, one learns that even 
rational functions are easier to deal with. Thus we start with one having a power 
of a quartic in its denominator. The evaluation of the identity 

f°° . . 

( - 1 ) J Q (x 4 + 2ax 2 + 1)"*+! = 2™+ 3 / 2 (a + 1)™+V2 m(a) ' 

where 

m 

(2.2) P m {a) = ^di(m)a' 
with 

(2.3) dj(m) = 2- 2m ^2 fc 

fc=; 

was first established in [7] . 



/ 2m — 2k\ ( m + k\ (k 



m J \ I 
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A standard hypergeometric argument yields 

(2.4) P m (a)=Pt^(a), 
where 

(2.5) P^\a) = p-ir- k (2 + I ) + * ^ + P ) 2"*(- + if 

is the classical Jacobi polynomial; the parameters a and (3 are given by a = m + | 
and j3 = —to — i. A general description of these functions and their properties 
are given in [T]. The twist here is that they depend on to, which means most of 
the properties of P m had to be proven from scratch. For instance, P m satisfies the 
recurrence 

(2TO-3)(4TO-3)a p (4m-3)fl(fl + l) , 

4to(to — l)(a — 1) 2to(to — l)(a — 1) 

4to(c 2 - 1) + 1 - 2a 2 

+ 7T1 T\ p m{a). 

2m(a — 1) 

This cannot be obtained by replacing a = m + \ and (3 = —to — | in the standard 
recurrence for the Jacobi polynomials. The reader will find in [3] several differnt 
proofs of (|2~T1) . 

The polynomials P m (a) makes a surprising appearance in the expansion 

(2.6) \J a + V\T~c = \/a+\ 



^ (-l) k P fc -!(a)c fc 
^ k 2 k + 1 {a + l) h 



as described in [9]. The special case a — 1 appears in [52], page 191, exercise 21. 
Ramanujan [4j had a more general expression, but only for the case c — a 2 : 



b k {n)a k 

[U T V J- T Lt~ J = J. -f 11a -|- J 

k=2 

where, for k > 2, 



(2.7) (a+Vl + a 2 ) n = l + na + J2 

k=2 



(2.8) 6 fc (n) 



n 2 (n 2 - 2 2 )(n 2 - 4 2 ) • • • (n 2 - (k - 2) 2 ) if k is even, 
n(n 2 - l 2 )(n 2 - 3 2 ) • • • (n 2 - (k - 2) 2 ) if k is odd. 



This result appears in Berndt [4] as Corollary 2 to Entry 14 and is machine- 
checkable, as are many of the identities in this section. 

The coefficients di(m) in (|2.3[) have :many interesting properties: 

• They form a unimodal sequence: there exists an index < to* < to such that 
dj (to) increases up to j = to* and decreases from then on. See [Hj for a proof of the 
more general statement: If P{x) is a polynomial with nondecreasing, nonnegative 
coefficients, then the coefficient sequence of P(x + 1) is unimodal. 

• They form a log-concave sequence: define the operator £({cifc}) := {a 2 — ak-ia^+i} 
acting on sequences of positive real numbers. A sequence {ak} is called log-concave 
if its image under £ is again a sequence of positive numbers; i.e. a 2 — ak-iak+i > 0. 
Note that this condition is satisfied if and only if the sequence {bk := log(afc)} is 
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concave, hence the name. We refer the reader to [53] for a detailed introduction. 
The log-concavity of di (to) was established in [35] using Computer Algebra tech- 
niques: in particular, cylindrical algebraic decompositions as developed in [21] and 

m- 

• They produce interesting polynomials: in one finds the representation 

(2-9) di(m)= „ ^",„ 

with 

m m 

(2.10) A l<m = on(m) Y[(^k - 1) - /9,(m) JJ(4A; + !)• 

fe=i fe=i 

Here aj and Pi are polynomials in m of degrees I and £ — 1, respectively. For 
example, ai(m) = 2m + 1 and P\(m) — 1, so that the coefficient of the linear term 
of P m (a) is 



-i / m m \ 

j^j (2m + l)n(4fc-l)-n( 4fc + 1 ) 

V fc=l fc=l / 



(2.11) di(m)= , 

m! 2 

J. Little established in [JUj the remarkable fact that the polynomials a; (to) and 
01 (to) have all their roots on the vertical line Re m = — i . 

When we showed this to Henry, he simply remarked: the only thing you have to 
do now is to let I — > oo and get the Riemann hypothesis. The proof in [40j consists 
in a study of the recurrence 

(2.12) y l+1 (s) = 2s m (s) - (s 2 - (21 - l) 2 ) y^s), 

satisfied by a/((s — l)/2) and fli((s — l)/2). There is no Number Theory in the 
proof, so it is not likely to connect to the Riemann zeta function ((s), but one never 
knows. 

The arithmetical properties of A^ m are beginning to be elucidated. We have 
shown that their 2-adic valuation satisfies 

(2.13) v 2 (A l;m ) =v 2 {{m + l-l) 2 i) + I, 

where (a)k = a(a + l)(a + 2) • • • (a + k — 1) is the Pochhammer symbol. This 
expression allows for a combinatorial interpretation of the block structure of these 
valuations. See [5] for details. 

3. The incipient rational landen transformation 

The clean analytic expression in (|2.1[) is not expected to extend to rational 
functions of higher order. In our analysis we distinguish according to the domain of 
integration: the finite interval case, mapped by a bilinear transformation to [0,oo), 
and the whole line. In this section we consider the definite integral, 

f*i\ nth ^ f°° cx^ + dx^ + e 
(3.1) Ue(a,b;c,d,e) = / -rr- dx, 



o x e + ax 4 + bx 2 + 1 

as the simplest case on [0,oo). The case of the real line is considered below. The 
integrand is chosen to be even by necessity: none of the techniques in this section 
work for the odd case. We normalize two of the coefficients in the denominator 
in order to reduce the number of parameters. The standard approach for the 
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evaluation of (J3JJ) is to introduce the change of variables x = tan 9. This leads to 
an intractable trigonometric integral. 

A different result is obtained if one first symmetrizes the denominator: we say 
that a polynomial of degree d is reciprocal if Qd{l/x) — x Qd(x), that is, the 
sequence of its coefficients is a palindrome. Observe that if Q d is any polynomial 
of degree d, then 

(3.2) T 2d (x) = x d Q d (x)Q d (l/x) 
is a reciprocal polynomial of degree 2d. For example, if 

(3.3) Q 6 (x) = x 6 + ax 4 + bx 2 + 1. 
then 

T12 (x) = x 12 + {a + b)x w + (a + b + ab)x 8 + 

+ (2 + a 2 + b 2 )x 6 + {a + b + ab)x 4 + (a + b)x 2 + 1. 

The numerator and denominator in the integrand of (|3.1|) are now scaled by x 6 Qe(l/x) 
to produce a new integrand with reciprocal denominator: 

(3.4) u^rp&La,, 

Jo T u {x) 

where we write 

5 6 

S 10 ( x ) = Y,s 3 x 2j and T 12 (ar) =^^a; 2j '. 

3=0 3=0 

The change of variables x — tan 9 now yields 

,„ TJ f /2 5io(tang) cos w {9) 

(3 - 5) U6 = Jo T u (t«atf) cos-(^) d °- 

Now let w — cos 29 and use sin 2 9 = |(1 — w) and cos 2 9 — |(1 +w) to check that 
the numerator and denominator of the new integrand, 



(3.6) Sio(tan 9) cos 10 9 = ^ s j s[n ~ J {i ( os 



and 



3=0 



6 

Ti 2 (tan 9) cos 12 6» = ^ tj sin 2j 9 cos 12 - 2j 9 

3=0 

= 2- 6 '£t j (i-wy(i + wf-J, 

3=0 

are both polynomials in w. The mirror symmetry of T12, reflected in = t6-ji 
shows that the new denominator is an even polynomial in w. The symmetry of 
cosine about ir/2 shows that the terms with odd power of w have a vanishing 
integral. Thus, with tp — 29, and using the symmetry of the integrand to reduce 
the integral from [0, 7r] to [0,7r/2], we obtain 

y _ r^ 2 r 4 cos 4 tp + r 2 cos 2 tp + r 

J Q q 6 cos 6 tp + q 4 cos 4 tp + q 2 cos 2 tp + q 
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The parameters Vj , qj have explicit formulas in terms of the original parameters of 
Uq . This even rational function of cos ip can now be expressed in terms of cos 2ip 
to produce (letting 8 <— 2ip) 

a.2 cos 2 8 + cti cos 8 + ccq 
6 = J Q f3 3 cos 3 8 + /3 2 cos 2 9 + 0~i cos 8 + J Q ' 

The final change of variables y = tan | yields a new rational form of the integrand: 
fov\ TT f°° C!y 4 + d iy 2 + ei 

Keeping track of the parameters, we have established: 



Theorem 3.1. The integral 

(3.8) U e = I fi "~" ' "~, : dx 

y ' Jo x 6 + ax 4 + bx 2 + 1 

is invariant under the change of parameters 

ab + 5a + 5b + 9 



(3.9) at 

bi 



{a + b + 2) 4 / 3 ' 
a + b + 6 



(a + fe + 2) 2 /3' 
for the denominator parameters and 

c + d + e 



Cl 

di 
ei 

for those of the numerator. 



(a + b+2) 2 / 3 ' 

(& + 3)c+2d+(a + 3)e 

a+b+2 ' 

c + e 



(a + b + 2) 1 / 3 



Theorem 13.11 is the precise analogue of the elliptic Landen transformation (jl.ip 
for the case of a rational integrand. We call (|3.9[) a rational Landen transformation. 
This construction was first presented in [8]. 



3.1. Even rational Landen transformations. More generally, there is a simi- 
lar transformation of coefficients for any even rational integrand; details appear in 
[lOj . We call these even rational Landen Transformations. The obstruction in the 
general case comes from (|3.6p ; one does not get a polynomial in w = cos 28. 

The method of proof for even rational integrals can be summarized as follows. 

1) Start with an even rational integral: 

[°° even polynomial in x 

U 2p = / : r-n dx. 

J even polynomial m x 
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2) Symmetrize the denominator to produce 

f°° even polynomial in x 

2p J Q even reciprocal polynomial in x 
The degree of the denominator is doubled. 

3) Let x = tan 9. Then 

f™^ 2 polynomial in cos 29 
2p J even polynomial in cos 29 

4) Symmetry produced the vanishing of the integrands with an odd power of cos 9 
in the numerator. We obtain 

f*/ 2 even polynomial in cos 29 
U 2p = I ; r-r-. — dO. 



(i 



even polynomial in cos 29 



5) Let ip = 26 to produce 



["* even polynomial in cos - !/" 

U 2p = / : — dtp. 

J even polynomial m cos ip 
Using symmetry this becomes an integral over [0, 7r/2]. 



6) Let y = tan-0 and use cos?/' = l/yl + y 2 to obtain 



DC 



. even polynomial in y 

U2 P = / ; r-r-. dy. 

Jq even polynomial m y 

The degree of the denominator is half of what it was in Step 5. 

Keeping track of the degrees one checks that the degree of the new rational func- 
tion is the same as the original one, with new coefficients that appear as functions 
of the old ones. 

4. A GEOMETRIC INTERPRETATION 

We now present a geometric foundation of the general even rational Landen 
transformation (]3.9[) using the theory of Riemann surfaces. The text |51j provides 
an introduction to this theory, including definitions of objects we will refer to here. 
The sequence of transformations in section [3] can be achieved in one step by relating 
tan 29 to tan#. For historical reasons (this is what we did first) we present the 
details with cotangent instead of tangent. 

Consider the even rational integral 

(4.1) 1=1 R{x)dx = -\ R{x)dx. 

Jo 2 J_ 00 



1 



Introduce the new variable 

(4-2) y = R 2 (x) = : 2( . 

motivated by the identity cot 29 = R 2 (cot 9) . The function R 2 
two-to-one map. The sections of the inverse are 



(4.3) x = a±{y) = y± yV + 1. 
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Splitting the original integral as 

R{x)dx + I R(x)dx 

-oo JO 

and introducing x = cr+(y) in the first and x = cr^(y) in the second integral, yields 

/oo 
(R+(y) + R-(y)) dy 
-OO 

where 

R+(y) = + R(a-(y)) and J?_(y) = JL^ (R(cr + (y)) - R(<r-(v)) ■ 

Vv + 1 

A direct calculation shows that i?+ and i?_ are rational functions of degree at most 
that of R. 

The change of variables y = R2{x) converts the meromorphic differential tp = 
R(x) dx into 

R(* + (y))^ + R(<r-(y)) d -^ = ((R(, + ) + R^)) + yiR ^\- R ^Ady 
dy dy y Vy 2 + l J 

= (R+{y)+R-(y)) dy. 

The general situation is this: start with a finite ramified cover it : X — > Y of 
Riemann surfaces and a meromorphic differential iponX. Let U <ZY he a simply 
connected domain that contains no critical values of 7r, and let <7i, . . . , Uk ■ U — > X 
be the distinct sections of it. Define 

k 

(4.6) n*ip T =y2<7j<p- 

u * — ' 

i=i 

In [33j we show that this construction preserves analytic 1-forms, that is, if cp is an 
analytic 1-form in X then ir^ip is an analytic 1-form in Y. Furthermore, for any 
rectifiable curve 7 on Y, we have 

(4.7) I K*<p= I V- 

In the case of projective space, this leads to the following: 

Lemma 4.1. If ir : P 1 — > P 1 is analytic, and ip = R{z) dz with R a rational 
function, then ir*tp can be written as R\{z) dz with R\ a rational function of degree 
at most the degree of R. 

This is the generalization of the fact that the integrals in (|4.ip and (|4. 5|) are the 
same. 

5. A FURTHER GENERALIZATION 

The procedure described in Section [3] can be extended with the rational map 
R m , defined by the identity 

(5.1) cot m9 = R m (cot 6). 

Here m € N is arbitrary greater or equal than 2. We present some elementary 
properties of the rational function R m . 
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Proposition 5.1. The rational function R m satifies: 

1) For m £ N define 

Lm/2J . . L(™-1)/2J / N 

Pm{x ) : = £ { -iy Qx— « and Q m (*) := g (jy + J^^- 

Then R m := P m /Q m . 

2) The function R m is conjugate to f m {x) := x m via M(x) := that is, 
Rm = M- 1 o f m o M. 

3) The polynomials P m and Q m have simple real zeros given by 

Pk := cot ( I for < k < m— f , and qk '■= cot ( — ) for 1 < k < m — 1. 

\ 2m / \m J 

If we change the domain to the entire real line, we can, using the rational sub- 
stitutions R m (x) i — ► a;, produce a rational Landen transformation for an arbitrary 
integrable rational function R(x) — B{x) / A{x) for each integer value of m. The re- 
sult is a new list of coefficients, from which one produces a second rational function 
RW(x) = J(x)/H(x) with 

(5.2) rm^rm*,. 



-oo A{X) J-oo H(x) 

Iteration of this procedure yields a sequence x n , that has a limit Xqq with conver- 
gence of order m, that is, 

(5.3) ||x„ +1 - Xoo || < Cllx^x^H" 1 . 

We describe this procedure here in the form of an algorithm; proofs appear in [41j . 

Lemma T4.II applied to the map ir(x) — R m (x), viewed as ramified cover of P 1 , 
guarantees the existence of a such new rational function R^\ The question of ef- 
fective computation of the coefficients of J and H is discussed below. In particular, 
we show that all these calculations can be done symbolically. 

• Algorithm for Deriving Rational Landen Transformations 

Step 1. The initial data is a rational function R(x) :— B{x) / A(x). We assume 
that A and B are polynomials with real coefficients and A has no real zeros and 
write 

(5.4) A{x) := E a kx p ~ k and B(x) := ^ hx p ~ 



-2-k 

k=0 fc=0 



Step 2. Choose a positive integer m > 2. 
Step 3. Introduce the polynomial 

(5.5) H(x) := Res z (A(z),P m (z)-xQ m (z)) 
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and write it as 

v 

(5.6) H(x) := J2 e l xP ~ 



1=0 



The polynomial H is thus defined as the determinant of the Sylvester matrix which 
is formed of the polynomial coefficients. As such, the coefficients e; of H(x) them- 
selves are integer polynomials in the a^. Explicitly 



p 

Mi 



(5.7) ei = (-l) l a^Y[Q m {xj) x af * {R m ( Xl ), R m (x 2 ), ...,R m {x p )), 



where {x\,x 2 , ■ ■ • ,x p } are the roots of A, each written according to multiplicity. 
The 

by 



(r>) 

The functions er\ are the elementary symmetric functions in p variables defined 



(5-8) ]J(y - yi) = Y j {-l) l a < f\y 1 , y p )y^ 1 . 

1=1 1=0 

It is possible to compute the coefficients e/ symbolically from the coefficients of A, 
without the knowledge of the roots of A. 
Also define 

(5.9) E{x) := H(R m {x)) x Q m {xf . 

Step 4. The polynomial A divides E and we denote the quotient by Z . The coef- 
ficients of Z are integer polynomials in the a^. 

Step 5. Define the polynomial C(x) := B(x)Z{x). 

Step 6. There exists a polynomial J(x), whose coefficients have an explicit formula 
in terms of the coefficients Cj of C(x), such that 

(5.10) r^-r#u 



A(x) ' J^H(x) 

This new integrand is the rational function whose existence is guaranteed by Lemma 
14.11 The explicit computation of the coefficients of J can be found in [U]. This is 
the rational Landen transformation of order m. 

Example 5.1. Completing the algorithm with m = 3 and the rational function 

1 



(5.11) R(x) = - 2 

ax z + ox + c 

produces the result stated below. Notice that the values of the iterates are ratios of 
integer polynomials of degree 3, as was stated above. The details of this example 
appear in [42] . 

Theorem 5.2. The integral 



12 



DANTE V. MANNA AND VICTOR H. MOLL 



is invariant under the transformation 

(5.13) a i — ► -| ((a + 3c) 2 - 36 2 ) , 

b - ^(3(a-c) 2 -6 2 ), 
c i ► ^ ((3a + c) 2 - 36 2 ) , 

where A := (3a + c)(a + 3c) — b 2 . The condition b 2 — 4ac < 0, imposed to ensure 
convergence of the integral, is preserved by the iteration. 



Example 5.2. In this example we follow the steps described above in order to 
produce a rational Landen transfromation of order 2 for the integral 



. , , , , b x 4 + hx 3 + b 2 x 2 + b 3 x + 6 4 

(5.14) 1 = / - = dx. 

_ oc aox° + a\X° + a^x^ + a^x A + a^x A + a^x + a§ 



Recall that the algorithm starts with a rational function R(x) and produces a 
new function £ 2 (i?(x)) satisfying 

/OO /'OO 
R(x)dx = l £ 2 (R(x))dx. 
-oo J — OO 



Step 1. The initial data is R{x) — B(x)/A(x) with 

(5.16) A(x) — aox e + aix 5 + a 2 x 4 + a 3 x 3 + a 4 x 2 + a 5 x + a 6 , 
and 

(5.17) B(x) = b x A + hx 3 + b 2 x 2 + b 3 x + b 4 . 
The parameter p is the degree of A, so p = 6. 

Step 2. We choose m — 2 to produce a method of order 2. The algorithm employs 
the polynomials P 2 {z) = z 2 — 1 and Q 2 (z) — 2z. 

Step 3. The polynomial 

(5.18) H(x) := Rcs z (A(z),z 2 - 1 - 2xz) 

is computed with the Mathematica command Resultant to obtain 

(5.19) H (x) = epx 6 + eix 5 + e 2 x 4 + e 3 x 3 + e 4 x 2 + e 5 x + e 6 , 



LANDEN SURVEY 



13 



with 

(5.20) 

e = 64a a 6 , 

ei = -32(a a 5 - ai<2 6 ), 

e2 = 16(ao<24 — aias + 6aoa6 + 0,20,6), 

e 3 = — 8(ao a 3 — 010,4 + 5oqo^ + a 2 as — baiae — 030,5), 

64 = 4(a a 2 - aia 3 + 4a a4 + 0204 — 4aia5 — a 3 a 5 + 9a a 6 + 4a 2 a 6 + ft4a 6 ), 

65 = — 2(aoai — a\02 + 3etoa 3 + 0203 — 3aio,4 — 0304 + 5oq05 + 

+3a2£i5 + 0405 — 5aiag — 3a 3 o,6 — o^oq), 

e 6 = ( a Q — a\ + a 2 — a 3 + a 4 — a 5 + a 6 )(a + a\ + a 2 + a 3 + a 4 + 05 + a 6 ). 

The polynomial H(x) is the denominator of the integrand £ 2 (i?(ir)) in (|5.15|) . 
In Step 3 we also define 

(5.21) E{x) = H(R 2 (x))Q 6 2 (x) = H (^^) ' C^xf . 

The function E(x) is a polynomial of degree 12, written as 



(5.22) E(x) = Y, 

a k x 



2-k 

l k x 

k=0 



Using the expressions for ej in (|5.20[) in (|5.2ip yields 

(5.23) ctQ = ol\2 = QAoqOq 

ai = — an — — 64(aoa5 — aide), 

oi2 = Oiio = 64(a a 4 — aia 5 + a 2 a 6 ), 

a 3 = -ag = -64(a a 3 - 0104 + a 2 a 5 - a 3 a 6 ), 

Q'4 = a 8 = 64(a a 2 — aia 3 + a 2 a 4 — a 3 a 5 + a 4 a 6 ), 

as = — 07 = — 64(aoai — a!a 2 + a 2 a 3 — a 3 a4 + 0405 — o^oq) 

a Ai 2 2,2 2,2 2 , 2\ 

a 6 = o4(a — a-y + a 2 — a 3 + a 4 — a 5 + a 6 ). 



Step 4. The polynomial A(a;) always divides E(x). The quotient is denoted by 
Z(a;). The values of a,- given in (|5.23p produce 

(5.24) Z(x) = 64(a — a\X + a 2 x 2 — a 3 x 3 + 04a; 4 — a 5 x 5 + OqX 6 ). 



Step 5. Define the polynomial C(x) := B{x)Z{x). In this case, C is of degree 10, 
written as 

10 

(5.25) C(x)=Y, c * xl0 ~ k > 

fe=0 
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and the coefficients Ck are given by 
(5.26) 



Co 


= 64a 6 6 , 






Cl 


= -64(a 5 6 


- a 6 bi), 




C2 


= 64(0460 — 


a$bi + o 6 6 2 ), 




C3 


= -64(a 3 6 


- a 4 6i + a 5 b 2 


- a 6 b 3 ), 


C 4 


= 64(a 2 feo - 


a 3 b\ + a 4 fe 2 - 


a$b 3 + a 6 & 4 ), 


C5 


= -64(ai& 


- a 2 bi + a 3 b 2 


- a 4 b 3 + a 5 6 4 ), 


C6 


= 64(a fo - 


aibx + a 2 b 2 - 


0363 + a 4 fe 4 ), 


C7 


= 64(a &i- 


a\b 2 + a 2 b 3 - 


a 3 6 4 ), 


C8 


= 64(a b 2 - 


aib 3 + a 2 bi), 




C9 


= 64(a 6 3 ~ 


ai& 4 ), 




ClO 


= 64ao&4. 







Step 6 produces the numerator J(x) of the new integrand £, 2 (R(x)) from the 
coefficients Cj given in (|5.26|) . The function J(x) is a polynomial of degree 4, 
written as 

4 

(5.27) J(x)=J2jkX i ~ k . 

k=0 

Using the values of (|5.26p we obtain 
(5.28) 

jo = 32(a 6 6 + a b 4 ), 

ji = -I6(a 5 b - a 6 bi + a b 3 ~ aib 4 ), 

j 2 = 8(a 4 & + 3a 6 6 - a 5 bi + aob 2 + a 6 b 2 - a\b 3 + 3a 6 4 + a 2 6 4 ), 
j 3 = -4(a 3 6 + 2a 5 6 + a-oh - a 4 6i - 2a 6 6i - a x b 2 + a 5 b 2 + 

+2a b 3 + a 3 b 3 ~ a 6 b 3 ~ 2a 1 b i - a 3 b 4 ), 
ji = 2(a b + a 2 b + a 4 b + a 6 b - a\bi - a 3 b\ - a 5 bi + a a b 2 + 

+a 2 b 2 + a 4 b 2 + a 6 b 2 - a x b 3 - a 3 b 3 — a 5 b 3 + a b 4 + a 2 b 4 + a 4 & 4 + a 6 6 4 ). 

The explicit computation of the coefficients of J can be found in [41] . 
The new rational function is 

(5.29) &2(R(x)):=^, 

with J(x) given in (|5.27p and H(x) in (|5.19p . The transformation is 

box 4 + hx 3 + b 2 x 2 + b 3 x + b 4 j x 4 + jix 3 + j 2 x 2 + j 3 x + j 4 



aox e + aix 5 + a 2 x 4 + a 3 x 3 + a 4 x 2 + a$X + egx 6 + e\x 5 + e 2 x 4 + e 3 x 3 + e 4 x 2 + 65a; + eg 

The numerator coefficients are given in (|5.20|) and the denominator ones in (|5.28|) . 
explicitly as polynomials in the coefficients of the original rational function. The 
generation of these polynomials is a completely symbolic procedure. 



LANDEN SURVEY 



15 



The first two steps of this algorithm, applied to the definite integral 

(5.30) [ 6 dX 3 = ^ (2\/3cos(7r/9) + Vz cos(27r/9) + 3 sin(27r/9)) 

produces the identities 

r°° 9f1 fii- 4 4- 1 9r 2 -4- 9r 4- 9) 

dx 



dx _ f 00 2(l6x' i + 12x 2 + 2x + 2) 



1 64a; 6 + 96a; 4 + 36a; 2 + 3 

4(2816a; 4 - 1024a; 3 + 8400a; 2 - 884a; + 5970) , 

dx 

12288a; 6 + 59904a; 4 + 87216a; 2 + 39601 

The convergence of the iterations of rational Landen transformations is discussed 
in the next section. 

6. The issue of convergence 

The convergence of the double sequence (a„, b n ) appearing in the elliptic Landen 
transformation (|1.1[) is easily established. Assume < bo < clq, then the arithmetic- 
geometric inequality yields b n < 6 n +i < a n+ i < a n . Also 

(6.1) < a„+i - b n +i = 77 ; — ; — 7=^- • 



2 0^+v^) 2 

This shows that a n and 6„ have a common limit: M = AGM(a, b), the arithmetic- 
geometric of a and b. The convergence is quadratic: 

(6.2) \a n+ i ~ M\ < C\a n - M\ 2 , 

for some constant C > independent of n. Details can be found in [15] . 

The Landen transformations produce maps on the space of coefficients of the 
integrand. In this section, we discuss the convergence of the rational Landen trans- 
formations. This discussion is divided in two cases: 

Case 1: the half-line. Let R(x) be an even rational function, written as R(x) = 
P(x)/Q(x), with 

p-i p 



a k x 2 ^ k \ 



(6.3) P(x) = J2 hx 2{p - 1 - k) , Q(x) = J2 

k=0 k=0 

and do = a p = 1. The parameter space is 

(6.4) «p+ = {(oi, • • • , Op_i; b , ■■■ , b p ^)} c W' 1 > 
We write 

(6.5) a := (ai, • • • , a p _i) and b := (b 0l • ■ • , 6 P ). 
Define 



(6.6) 



/>oo 

A2 P = {(ffli, • • • , a p _i) £ : / dx is finite}. 

Jo 



Observe that the convergence of the integral depends only on the parameters in the 
denominator. 

The Landen transformations provide a map 

(6-7) $ 2p : ?P+ - 5JJ+ 
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that preserves the integral. Introduce the notation 

(6.8) a n = (a<"\ • • • , a { fy and b„ = (&£ n) , ■ • • 
where 

(6.9) (a„,b„) = $ 2 j,(a„_i,b„_i) 
are the iterates of the map $2 P - 

The result that one expects is this: 
Theorem 6.1. The region A 2p is invariant under the map $ 2 p- Moreover 



(6.10) 



and there exists a number L, that depends on the initial conditions, such that 

( 6 .iD (CTM'T >•■■■(;:!»■ 

This is equivalent to say that the sequence of rational functions formed by the 
Landen transformations, converge to L/(x 2 + 1). 

This was established in |33j using the geometric language of Landen transforma- 
tions which, while unexpected, is satisfactory. 

Theorem 6.2. Let ip be a 1-form, holomorphic in a neighborhood oft C P 1 . Then 
(6.12) " " ' ' ' ' "' ; 



lim (ir.) n (p = - ( f tp) — 



where the convergence is uniform on compact subsets of U, the neighborhood in 
the definition of 7r* . 

The proof is detailed for the map n(z) = z 2 ~ 1 = R,2(z), but it extends without 
difficulty to the generalization R m . 

Theorem 16.21 can be equivalently reformulated as: 

Theorem 6.3. The iterates of the Landen transformation starting at (ao,bo) £ 
converge (to the limit stated in Theorem 16. 11 ) if and only if the integral formed 
by the initial data is finite. 

It would be desirable to establish this result by purely dynamical techniques. 
This has been established only for the case p = 3. In that case the Landen trans- 
formation for 

cx 4 + dx 2 + e 
x 6 + ax 4 + bx 2 + 1 



IS 



,„.,,, ab + 5a + 5b + 9 

(6 ' 14) ^ ~ ( fl + b + 2) 4 /3 ' 

a + b + 6 

bi 



(a + fe + 2) 2 /3' 
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Figure 1. Zero locus of R(a,b) 



coupled with 



The region 
(6.15) 



Cl 

di 
ei 



(a + fe+2) 2 / 3 ' 

(b + 3)c + 2d+ (a + 3)e 



a + b + 2 

c + e 



(a + fe+2) 1 /3- 



A 6 = {(a, 6) £l 2 :[/ 6 <oo} 



is described by the discriminant curve the zero set of the polynomial 



(6.16) 



R(a, b) = Aa A + Ab A - 18ab - a 2 b A + 27 



This zero set, shown in Figure [TJ has two connected components: the first one 9t + 
contains (3, 3) as a cusp and the second one given by i?_(a, b) = 0, is disjoint 
from the first quadrant. The branch 9l_ is the boundary of the set Aq. 
The identity 

(a - bf R(a, b) 



(6.17) 



R(a u h) 



(a + 6 + 2) 4 ' 

shows that is invariant under $ 6 . By examining the effect of this map along 
lines of slope —1, we obtain a direct parametrization of the flow on the discriminant 
curve. Indeed, this curve is parametrized by 



(6.18) 



a{s) = 



s A + 4 



and b(s) 



s J + 16 
4s 
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Then 

'4(s 2 + 4) 



2 i A\2\ 1/3 



(6-19) ¥>(•)-, fi(s + 2)2 

gives the image of the Landen transformation <i>6; that is, 

(6.20) $ 6 (a(«), 6(a)) = (afcOO), %(«)). 

The map $6 has three fixed points: (3, 3), that is super-attracting, a saddle point 
Pi on the lower branch 9t_ of the discriminant curve, and a third unstable spiral 
below this lower branch. In [25] we prove: 

Theorem 6.4. The lower branch of the discriminant curve is the curve Aq. This 
curve is also the global unstable manifold of the saddle point Pi. Therefore the 
iterations of $6 starting at (a, b) converge if and only if the integral Us, formed 
with the parameters (a, b), is finite. Moreover, (a n ,b n ) — > (3,3) quadratically and 
there exists a number L such that (c n , d n , e n ) —* (1, 2, 1)L. 

The next result provides an analogue of the AGM (|1.12p for the rational case. 
The main differences here are that our iterates converge to an algebraic number 
and we achieve order-m convergence. 

Case 2: The whole-line: This works for any choice of positive integer m. Let 
R(x) be a rational function, written as R(x) = B(x)/A(x). Assume that the 
coefficients of A and B are real, that A has no real zeros and that deg(B) < 
deg(A) — 2. These conditions are imposed to guarantee the existence of 

P oo 

(6.21) 1=1 R(x)dx. 



(6.22) A(x) = a k x p ~ k and B{x) = ^ b k x* 



In particular A must be of even degree, and we write 

P P-2 

r p-2-k 

'K- 

k=0 k=Q 

We can also require that deg(gcd(A, B)) = 0. 

The class of such rational functions will be denoted by 9\ p . 

The algorithm presented in Section [5j provides a transformation on the parame- 
ters 

(6.23) ¥ P := {a , a u --- ,a p ;b , b u --- , 6 p _ 2 } = x M^ 1 

of R S 9\ p that preserves the integral I. In fact, we produce a family of maps, 
indexed by to G N, 

^m,p ■ ^p y 

such that 

/oo />oc 
R[x)dx = I £ m<p (R(x)) dx. 
-oo J — OO 

The maps £ m , p induce a rational Landen transformation on the parameter space: 

(6.25) $ m ,p : <Pp - <Pp 

by simply listing the coefficients of the function £, m ^ p (R(x)). 
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The original integral is written in the form 

(6.26) i = h. r xp - 2 + ft °; 6 ^- 3 + & ° ;^- 4 + ■ ■ ■ + »o-y » dx . 

a o J~oo x p + % aixP 1 + a a 2 xP 2 + ■ ■ ■ + a a p 
The Landen transformation generates a sequence of coefficients, 

(6.27) qj p>n := {a^UfV^aW;^,^,.-',^}. 
with typfi — ^P P as in l|6.23[) . We expect that, as n — > oo, 



(6.28) 
converges to 



J °0 °0 & & & 



(6.29) Xoo := I 0, (*) ,0,( q l---,( 9 \0,( q : 1 ),0/ q -^ ^ 



T/ V2y y^y V 1 y V 2 y ~ 1 

where q = p/2. Moreover, we should have 
(6.30) Hx^i - xooll < C||x„- Xoo |r. 

The invariance of the integral then shows that 

(6-31) - -I. 

This produces an iterative method to evaluate the integral of a rational function. 
The method's convergence is of order-m. 

The convergence of these iterations, and in particular the bound (|6.30|) . can be 
established by the argument presented in sectionHJ Thus, the transformation Z m ,p 
leads to a sequence that has order-m convergence. We expect to develop these ideas 
into an efficient numerical method for integration. 

We choose to measure the convergence of a sequence of vectors to is in the 
L 2 — norm, 



V2 

( w„,\\A 

^2p~2 



2p-2 N 
2 



(6.32) Wa = ^(Ell^ , 

. fc=i / 

and also the Loo-norm, 

(6.33) IMU = Max{|| Vfe || :\<k<2p-2}. 

The rational functions appearing as integrands have rational coefficients, so, as a 
measure of their complexity, we take the largest number of digits of these coeffi- 
cients. This appears in the column marked size. 

The following tables illustrate the iterates of rational Landen transformations of 
order 2, 3 and 4, applied to the example 

Ux 3 + 74a; 2 + 184.x + 208 °° ~ ~ 12' 



2(1 



DANTE V. MANNA AND VICTOR H. MOLL 



The first column gives the L2-norm of u n — Woe, the second its Loo-norm, the third 
presents the relative error in (|6.31j) . and in the last column we give the size of the 
rational integrand. At each step, we verify that the new rational function integrates 
to -7tt/12. 

Method of order 2 



n 


L 2 -norm 


-norm 


Error 


Size 


1 


58.7171 


69.1000 


1.02060 


5 


2 


7.444927 


9.64324 


1.04473 


10 


3 


4.04691 


5.36256 


0.945481 


18 


4 


1.81592 


2.41858 


1.15092 


41 


5 


0.360422 


0.411437 


0.262511 


82 


6 


0.0298892 


0.0249128 


0.0189903 


164 


7 


0.000256824 


0.000299728 


0.0000362352 


327 


8 


1.92454 x 10~ 8 


2.24568 x 10~ 8 


1.47053 x 10~ 8 


659 


9 


1.0823 x 10~ 16 


1.2609 x 10~ 16 


8.2207 x 10~ 17 


1318 



As expected, for the method of order 2, we observe quadratic convergence in the 
L2— norm and also in the Loo— norm. The size of the coefficients of the integrand 
is approximately doubled at each iteration. 



Method of order 3 



n 


L2-norm 


Loo-norm 


Error 


Size 


1 


15.2207 




20.2945 


1.03511 




8 


2 


1.97988 




1.83067 


0.859941 




23 


3 


0.41100 




0.338358 


0.197044 




69 


4 


0.00842346 




0.00815475 


0.00597363 




208 


5 


5.05016 x 10" 


-8 


5.75969 x 10~ 8 


1.64059 x 10" 


-9 


626 


6 


1.09651 x 10- 


23 


1.02510 x 10~ 23 


3.86286 x 10- 


24 


1878 


7 


1.12238 x 10- 


70 


1.22843 x 10~ 70 


8.59237 x 10- 


71 


5634 


Method of order 4 


71 


L2-norm 


Loo -norm 


Error 


Size 


1 


7.44927 




9.64324 


1.04473 




10 


2 


1.81592 




2.41858 


1.15092 




41 


3 


0.0298892 




0.0249128 


0.0189903 




164 


4 


1.92454 x 10" 


-8 


2.249128 x 10~ 8 


1.47053 x 10- 


-8 


659 


5 


3.40769 x 10" 


33 


3.96407 x 10~ 33 


2.56817 x 10- 


33 


2637 



Example 6.1. A method of order 3 for the evaluation of the quadratic integral 

/OO j 
_oo A£ + OX + C 

has been analyzed in [32] • We refer to Example 15.11 for the explicit formulas of this 
Landen transformation, and define the iterates accordingly. From there, we prove 
that the error term, 

(6.35) e„ := (a n - \\J Aac - b 2 , b n , c n - \ \/ Aac - b 2 j 
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satisfies e„ — > as n — > oo, with cubic rate: 

(6.36) ||e n+ i|| <q|e„|| 3 . 

The proof of convergence is elementary. Therefore we have 



(6.37) («„, b n , c„) -> (Vac-6 2 /4, 0, ^4ac - 6 2 /4) , 
which, in conjunction with (|6.34p , implies that 

2 f°° dx 

(6.38) I 



V4ac — b 2 J-oo x 2 + 1 

exactly as one would have concluded by completing the square. Unlike completing 
the square, our method extends to a general rational integral over the real line. 



7. The appearance of the AGM in diverse contexts 
The (elliptic) Landen transformation 

(7.1) a\ <— — (a + b) and b\ <— \fab~ 

leaving invariant the elliptic integral 

r^ 2 d v 



(7.2) G(a,b) = 



\J a 2 cos 2 ip + b 2 sin 2 ip 



appears in many different forms. In this last section we present a partial list of 
them. 



7.1. The Elliptic Landen Transformation. For the lattice L = Zffi loZ, intro- 
duce the theta-functions 

oo 

(7.3) $ 3 (x,u>):= z 2n q n2 &nddi(x,uj) 

n— — oo 

where z — e mx and q — e™. The condition Imcj > is imposed to ensure 
convergence of the series. These functions admit a variety of remarkable identities. 
In particular, the null-values (those with x = 0) satisfy 

tf 2 (0, 2w) = tf 3 (0, w)0 4 (O, w) and ?? 2 (0, 2uj) = 1 (^|(0, u) + #1(0, to)) , 

and completely characterize values of the AGM, leading to the earlier result [T5] , 
Grayson [3U] has used the doubling of the period u> to derive the arithmetic- 
geometric mean from the cubic equations describing the corresponding elliptic 
curves. See Chapter 3 in [43] for more information. P. Sole et al. [49l [50] have pro- 
duced proofs of generalizations of these identities by lattice enumeration methods 
related to binary and ternary codes. 



DC 
71 — — OO 
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7.2. A time-one map. We now present a deeper and more modern version of a 
result known to Gauss: given a sequence of points {x n } on a manifold X , decide 
whether there is a differential equation 

(7.4) g = V(x), 

starting at xo such that x n — x(n,xo). Here x(t,xo) is the unique solution to (|7.4p 
satisfying x(0,xq) = Xq. Denote by 

(7.5) (l> e iiip(a,b) = Q(a + &), Vab\ 

the familiar elliptic Landen transformation. Now take a,b GM. with a > b > 0. Use 
the null- values of the theta functions to find unique values (r, p) such that 

(7.6) a = pi?|(0,r) and6 = /9i?|(0,T). 
Finally define 

(7.7) x ellip (t) = (o(t), &(*)) = P (^(0, 2V), i?*(0, 2*t)) , 
with a; e iji p (0) = (a, 6). The remarkable result is [28] : 

Theorem 7.1. (Deift, Li, Previato, Tomei). The map £ — > x e iu p (t) is an integrable 
Hamiltonian flow on X equipped with an appropriate symplectic structure. The 
Hamiltonian is the complete elliptic integral G(a, b) and the angle is (essentially the 
logarithm of) the second period of the elliptic curve associated with r. Moreover 

(7-8) x e iup(k) = 4> k elKp (a,b). 

Thus the arithmetic-geometric algorithm is the time-one map of a completely inte- 
grable Hamiltonian flow. 

Notice that this theorem shows that the result in question respects some addi- 
tional structures whose invention postdates Gauss. 

A natural question is whether the map (|3.Q|) appears as a time-one map of an 
interesting flow. 



7.3. A quadruple sequence. Several variations of the elliptic Landen appear 
in the literature. Borchardt [5] considers the four-term quadratically convergent 
iteration 

, n n x a n + b n + c n + d n 

(7.9) a n+ i 

bn+l 
Cn+1 





4 


V a n b n 


+ V c n d n 




2 




+ \Jb n d n 




2 


V O-ndn 


+ \Jb n c n 



2 

starting with qq = a, 6q — b, cq — c and o?o = d. The common limit, denoted by 



d n +i 

i 

G{a, b, c, d), is given by 

(710) 1 _ 1 [" '■ r I'' -!!)<>■'■ <hl 



G{a,b, Cl d) tt 2 7 J ai y /R( x )R(y) ' 
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where R(x) — x(x — ao)(x — ai)(x — 0L2){x — as) and the numbers ay are given by 
explicit formulas in terms of the parameters a, b, c, d. Details are given by Mestre 
in 05. 

The initial conditions (a, 6, c, d) £ R 4 for which the iteration converges has some 
interesting invariant subsets. When a = b and c = d, we recover the AGM iteration 
(|1.1|) . In the case that b = c = d, we have another invariant subset, linking to an 
iterative mean described below. 



7.4. Variations of AGM with hypergeometric limit. Let A € N. The anal- 
ysis of 

. a„ + (AT - 1)6„ a„ - 6„ 

(7.11) a n+ i = — and c„+i 



A " T± A ' 

with b n = (a„ — c^) 1 ^ , is presented in [TB]. All the common ingredients appear 
there: a common limit, fast convergence, theta functions and sophisticated itera- 
tions for the evaluation of ir. The common limit is denoted by AGAr(a,6). The 
convergence is of order A and the limit is identified for small A: for < k < 1, 



1 



AG 2 (l,fc) 
1 

AG 3 (l,fc) 



2 Fi(l/2, 1/2; 1; 1 - k 2 ) and 



2^(1/3, 2/3; 1; 1 - k 2 ) 



where 

00 

(7.12) 2 F 1 (a, b; c; x) = ^ 



(a) k ( b )k x k 



is the classical hypergeometric function. There are integral representations of these 
as well which parallel (|1.12p . see [13], Section 6.1 for details. 

Other hypergeometric values appear from similar iterations. For example, 

(7.13) a n+ i = a " + and b ^ _ ^/ bn ( an + b n )/2, 
have a common limit, denoted by A^a, b). It is given by 

(7.14) ^ = ^(1/4,3/4; l;l-k 2 ). 

To compute tt quartically, start at a = 1, 6 = (12\/2 — 16) 1 / 4 . Now compute 
a n from two steps of AG 2 : 



(7.15) a n+ i = , and b n+i 

Then 



(7.16) 7r = Jim 3a 4 +1 ( 1 - £ 2'+ 1 (a 4 - a 4 +1 ) 
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with |<2n+l — 7r| < C\a n 
the partial sums of 

(7.17) tt = 4^ ( l) 



^ 2k 

k=0 



The sequence (a n , b n ) defined by the iteration 

1 

, 7rs , a n + 2fe ra ( b n {al + a n b n + b 2 n ) \3 
l'- 18 J a n+ i = , b n+ i = I I , 

starting at ao = 1, 60 = x are analyzed in [17) . They have a common limit -F (a;) 
given by 

(7.19) _L_ =2Fl (l,|;l;l. 



7.5. Iterations where the limit is harder to find. J. Borwein and P. Borwein 
[16] studied the iteration of 

,_ ofV . , ,s (a + 36 ^06 + 6^ 
(7-20) (a,6)->l^— , I, 

and showed the existence of a common limit B(ao, bo). Define B(x) — B(l,x). The 
study of the iteration (|7.20[) is based on the functional equation 

(7.21) B( X ) = 

and a parametrization of the iterates by theta functions |16| . The complete analysis 
of (|7.20p starts with the purely computational observation that 

(7.22) B(x) ~ — log _2 (x/4) as x -> 0. 

o 

H. H. Chan, K. Chua and P. Sole [21] identified the limiting function as 
(7-23) B{ X )=(,F 1 (l,\-l;27f^$)y\ 

valid for | < x < 1. A similar hypergeometric expression gives B(x) for < x < I. 



7.6. Fast computation of elementary functions. The fast convergence of the 
elliptic Landen recurrence (|1.1|) to the arithmetic-geometric mean provides a method 
for numerical evaluation of the elliptic integral G(a, b). The same idea provides for 
the fast computation of elementary functions. For example, in [14] we find the 
estimate 

(7.24) \\ogx - (G(l, 1(T") - G(l, KT"?)) | < tiKT 2 ^ 1 ), 

for < x < 1 and n > 3. 



LANDEN SURVEY 



25 



7.7. A continued fraction. The continued fraction 

(7.25) R n (a, b) = a — 2 

Tj H 



4a 2 



9b 2 



!? + ••• 

has an interesting connection to the AGM. In their study of the convergence of 
R n (a,b), J. Borwein, R. Crandall and G. Fee [TS] established the identity 

(7.26) R v VaM^j = i (R v (a, b) + R n (b, a)) . 

This identity originates with Ramanujan; the similarity with AGM is now direct. 

The continued fraction converges for positive real parameters, but for a, b G C 
the convergence question is quite delicate. For example, the even/odd parts of 
J?i(l,i) converge to distinct limits. See [19] and [20] for more details. 



7.8. Elliptic Landen with complex initial conditions. The iteration of (jl.ip 
with ao, bo G C requires a choice of square root at each step. Let a, b £ C be 
non-zero and assume a ^ ±6. A square root c of ab is called the right choice if 



(7.27) 



< 



a + b 



+ c 



It turns out that in order to have a limit for (jl.ip one has to make the right choice 
for all but finitely many indices n > 1. This is described in detail by Cox [27] . 

7.9. Elliptic Landen with p-adic initial conditions. Let p be a prime and a, b 
be non-zero p-adic numbers. In order to guarantee that the p-adic series 



(7.28) 



i=0 



-1 



converges, and thus defines a p-adic square root of ab, one must assume 
(7.29) b/a=lmodp a , 

where a = 3 for p = 2 and 1 otherwise. The corresponding sequence defined by 
(II. ip converges for p ^ 2 to a common limit: the p-adic AGM. In the case p = 2 
one must assume that the initial conditions satisfy b/a = 1 mod 16. In the case 
b/a = 1 mod 8 but not 1 modulo 16, the corresponding sequence (a„, b n ) does not 
converege, but the sequence of so-called absolute invariants 

. VK-albl + bjf 
{ ] Jn ~ <bl(al~blY 

converges to a 2-adic integer. Information about these issues can be found in [32] . 
D. Kohel [3S] has proposed a generalization of the AGM for elliptic curves over 
a field of characteristic p e {2, 3, 5, 7, 13}. Mestre [IS] has developed an AGM 
theory for ordinary hyperelliptic curves over a field of characteristic 2. This has 
been extended to non- hyperelliptic curves of genus 3 curves by C. Ritzenhaler [38] . 
An algorithm for counting points for ordinary elliptic curves over finite fields of 
characteristic p > 2 based on the AGM is presented in R. Carls [23] . 
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7.10. Higher genus AGM. An algorithm analog to the AGM for abelian integrals 
of genus 2 was discussed by Richelot [47], [48] and Humberdt [34] . Some details 
are discussed by J. Bost and J. F. Mestre in [5T]. The case of abelian integrals of 
genus 3 is due to D. Lehavi [38] and D. Lehavi and C. Ritzenthaler 39J. 

Gauss was correct: his numerical calculation (|1.4p has grown in many unexpected 
directions. 
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